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The two dimensional XY spin glass in is studied numerically by a finite size defect energy scaling 
method at T = in the vortex representation which allows us to compute the exact (in principle) 
spin and chiral domain wall energies. We confirm earlier predictions that there is no glass phase at 
any finite T. Our results strongly support the conjecture that both spin and chiral order have the 
same correlation length exponent i^s = Vc ~ 2.70. Preliminary results in 3d are also obtained. 
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The XY spin glass has been the subject of consider- 
able attention and controversy for some time and is still 
not understood. It has been known since the seminal 
work of Villain ||] that vector spin glass models have 
chiral or reflection symmetry in addition to the contin- 
uous rotational symmetry. Consequently, the XY spin 
glass may have two different glass orders, a spin glass 
order and a chiral glass order. It is widely accepted and 
has become part of the spin glass folklore that, in two 
and three dimensions, chiral and spin variables decouple 
at long distances and order independently (|-|^ although 
there is a hint that this may not hold in four dimensions 
Numerical estimates of the correlation length expo- 
nents V in two dimensions, where both spin and chiral 
order set in at T = as ^^.c ^ T~'^= '=, indicate that 
i/c = 2.57 ± 0.003 and = 1.29 ± 0.02 fe) which agree 
with older, less accurate estimates ||,^,^]- The decou- 
pling of chiral and spin degrees of freedom seems to be 
well established by these numerical results, but some an- 
alytic work on special models ||8|~|l0t] implies that, at least 
for XY spin glasses below their lower critical dimension 
di > 2 when order sets in at T = 0, both correlation 
lengths diverge with the same exponent Vs = i^c- To add 
to the confusion, there is rather convincing evidence that 
chiral order sets in for < T < Tc, while spin glass order 
occurs only at T = in Sd [p|-^|JTl|-[r^ . These numeri- 
cal investigations have led to another piece of accepted 
folklore, namely that the lower critical dimension di > 4 
for spin glass order [|l3|,^. A very recent simulation H] 
concluded that earlier simulations are misleading because 
the spin defect energy began to grow with system size L 
at values of L just beyond the limit accessible to earlier 
attempts and that di is slightly less than three. However, 
chiral order is robust in 3d. In 2d, all simulations agree 
that chiral and spin glass order set in at T = but with 
different exponents Vc ~ "^i^s ~ 2.6. 

The theoretical situation is unclear since, to our knowl- 
edge, there is no unambiguous proof of any of the ac- 
cepted folklore outlined above ||l^-[l^, numerical simu- 
lations are contradictory j^] and the analytic work on 
special models |8|-|l0| is difficult to reconcile with the ap- 
parently unambiguous numerical simulations on the 2d 



XY spin glass. In this letter, we attempt to clarify the 
contradictory conclusions from numerical and analytic 
studies outlined above and to identify which should be 
retained and which need revision. Our essential conclu- 
sion is that, by carefully defining spin and a chiral domain 
wall energies, we find numerical agreement with the con- 
jecture that 9s = 6c in 2d where 0s,c — —^/vs^c are 
the T = stiffness exponents. Although the conjecture 
is not rigorous, it is the only, to our knowledge, analytic 
prediction existing and is the only check we have on the 
validity or otherwise of the numerical method used, at 
least until some rigorous testable predictions are made. 
If one accepts that a valid numerical simulation must 
agree with the conjecture, the implications go far beyond 
minor points such as the numerical values of stiffness ex- 
ponents but implies that most of the XY spin and gauge 
glass folklore is incorrect. The lower critical dimension 
for both spin glass and chiral order is 2 < < 3, the 
chiral glass scenario 9s < Q and 0c > in id is not pos- 
sible but both stiffness exponents are positive and the 
presently accepted numerical values in 2d and id are in- 
correct and need re-examination. 

A natural way of investigating order is to compute the 
domain wall or defect energy AE{L) of a system of size L 
for several realizations of disorder (samples) for different 
values of L and fit to the finite size scaling ansatz [p^jlS 



< l\E{L) >- 



(1) 



where < • ■ • > denotes an average over disorder, 
AE{L) — Eo{L) — Eq{L) the domain wall energy with 
Eo{L) the ground state (GS) energy, E£i{L) the energy 
of the system of size L containing a spin or chiral domain 
wall and Os,c is the spin (s) or chiral (c) stiffness expo- 
nent. There are two main difficulties in applying these 
ideas to a finite disordered system. The first is how to 
define Eq and E^ for a finite system with disorder since 
the GS configuration is unknown and the energy of a fi- 
nite system depends on the boundary conditions (BG) 
imposed which must be compatible with the GS config- 
uration. A spin or chiral domain wall is induced by an 
appropriate change in these BC and Ejj is the minimum 
energy of the system subject to these new BC. The sec- 
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ond difficulty is the computational problem of finding Eq 
and Ed sufficiently accurately so the errors in AE{L) 
can be controlled and kept small. The numerical data is 
fitted to eq.(§) in an attempt to verify the scaling ansatz 
and to obtain numerical values of the fundamental stiff- 
ness exponents 9s and 9c- These constraints limit the 
accessible sizes L to small values when the BC have large 
effects and it is essential to treat the BC properly to de- 
fine Eq and Ed consistently for a fit of the numerical 
data to eq. (|^) to have any meaning. 

The Hamiltonian of a ±J XY spin glass on a L x L 
square lattice is 

H=J2V{0,^9,-A,,) (2) 

<ij> 

where T^((/)) is an even 27r periodic function of (p with 
a maximum at (/) = tt, usually taken to be V(</'ij) = 
—JijCos{(j)ij) with the coupling Jij = J > for ij near- 
est neighbor sites of a square lattice. The random bond 
variables Aij = 0, tt with equal probability 1/2 corre- 
spond to ferro and antiferromagnetic coupling between 
neighboring spins. We imagine the system of eq.(|^) on a 
torus which corresponds to imposing periodic BC on the 
phases 9i^^i^ = ^i^+L^y = 9i^^i^+L with ia:^y = l,---,L 
and coupling spins on opposite faces by some interaction 
V{9L,iy,9i^i^) and V{9i^,L,9i^^i) which may be regarded 
as defining the BC. In principle, the GS is obtained by 
minimizing the energy with respect to the bulk vari- 
ables 9i and all possible V. This program is too difhcult 
for us so we will restrict ourselves to V which impose 
a spin defect or a chiral defect. To impose a spin do- 
main wall, we choose V = V{9i — 9j — Aij) where the 
phase differences Aij between corresponding sites i,j on 
opposite faces may be varied to find the minimum en- 
ergy Eq. It is not necessary to vary every Aij as each 
elementary plaquette on the torus is equivalent and the 
plaquettes between opposite faces are indistinguishable 
from the others and play no special role. We therefore 
keep fixed the frustrations /r = X^nr /^i" fixed where 
the sum is over the bonds in a clockwise direction of the 
elementary plaquette whose center is at the site r of the 
dual lattice. We are free to choose V to impose a global 
phase twist = 0, tt in the direction /i round the torus. 
The lowest energy i?o(A^) is 27r periodic in with a 
minimum at some AJJ which depends on the particular 
sample. To introduce a spin domain wall perpendicular 
to X, one simply changes the twists from their best twist 
(BT) values A° (A^ + tt, A°) and find the minimum 
energy subject to this constraint, which yields the energy 
with a spin domain wall Eajj^L) > Eo{L). The spin de- 
fect energy AEf'^{L) = Esd{L) — Eq{L) is computed for 
different samples and sizes L and fitted to 

< Ai?f^(L)>^L^?" (3) 

to obtain the spin stiffness exponent 6f^. A chiral do- 
main wall is imposed by reflective BC [0M| which means 



that there is a seam encircling the torus in (say) the y di- 
rection across which the spins interact as ^ = V{9i+9j — 
Aij ) which is equivalent to a reflection of the spins about 
some arbitrary axis. In principle, one can follow the same 
procedure as for the spin domain wall to obtain the chiral 
defect energy AEf'^{L) = Ecd{L) - Eo{L) where E^d 
is the minimum energy with the modified interactions 
on a seam. However, there is no reason to expect that 
EcD > Eq as the BC defining Eq may trap a chiral defect 
in some samples in which cases the modified interactions 
cancel the chiral defect and Ecd < Eq, as in fact does 
occur. We therefore define AE^'^ =| S^d -£^o I, average 
this over disorder and fit to < AEf'^ >^ L^" to obtain 
the chiral stiffness exponent. This does not affect AEf^ 
as both EsD and Eq contain the same chiral defects. The 
procedure described above using the phase representation 
of eq. (H) is similar to that of most previous studies 
except that these omit the minimization with respect to 
the twists A^, apply naive periodic and antiperiodic BC 
and call the lowest energies Ep and Eap- Neither of these 
BC is compatible with the GS and both must induce 
some excitation from Eq. Nevertheless, the spin defect 
energy is defined as AEf^ =| Eap — Ep \ and the spin 
stiffness exponent defined by < AEf'^{L) >~ . We 
call this a random twist (RT) measurement as both BC 
are equivalent to some random choice of A^ relative to 
A° for each sample. There is no good reason to expect 
AE^'^{L) to scale as L^" but if it does, there is less rea- 
son to expect any relation between 9f^^ and 9f^ or 9c. 

The procedure in terms of the phase representation of 
the XY spin glass Hamiltonian of eq.(||) is followed by 
previous studies. The aim is to obtain AE{L) by inde- 
pendently minimizing the Hamiltonian with respect to 
the 9i to obtain Eo and Eq. This requires finding es- 
sentially exact global minima for each sample to control 
the errors in < AE{L) > to be purely statistical and 
0(iV~^/^) where N is the number of samples. If the min- 
imization algorithm fails to find the true global minima, 
the errors in < AE{L) > will be uncontrolled and very 
large, making the data point useless. Since the \ 0i \< tt 
are continuous, one has to perform a numerical search 
of a huge configuration space, most of which does not 
even correspond to a local energy minimum. To reduce 
the volume of the space, we transform to a Coulomb gas 
(CG) representation which eliminates spin wave excita- 
tions and parametrizes the problem in terms of integer 
valued vortex or charge configurations, each of which is 
a local energy minimum. This reduces the space to be 
searched to manageable size although it introduces long 
ranged Coulomb interactions between vortices. The po- 
tential V{4>) in cq.(||) is taken as a piecewise parabolic 
potential equivalent to a Villain potential at T = 

<ij> 
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where = —nji is any integer on the bond ij. By a 
duality transformation 20,2^] the Coulomb gas Hamil- 
tonian with periodic BC in the phases becomes 

2^V^(<7r-/r)G(r-r')(gr'-/r') 



where 



+ J{al+al)/2L^ 
ax = -27r[L{qxi - f^i) + ^(^r - U)y\ 

r 

Gy = -2Tr[L{qyi - fyi) - ^{qr - fr)x] 

r 



(5) 



/-^ 4 — 2coskx — 2cosky 



(6) 
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Here, r = {x, y) denotes the sites of the dual lattice 
and G(r) is the lattice Green's function. In eq.(H), 
ka = 2'Kna/ L with = (0, 1, • • ■ , L — 1). The topo- 
logical charge, qr^ is the circulation of the phase about 
the plaquette at r and can be any integer subject to the 
neutrality condition q^. ~ 0. The frustration at r, 
/r = X^Dr IS the circulation of Aij round the pla- 

quette. fxi = X^Da; ^ii/^"" is the circulation round the 
whole torus on the x bonds of plaquettes &t y — 1 and q^i 
is the circulation of the phase, fyi and qyi are defined 
similarly. Periodic BC in the phases 6i restrict qxijqyi 
to be integers. A chiral domain wall is introduced by 
reflective BC when the Hamiltonian becomes [|| 



Hn = 27r2 jE(<Zr - /r)(gr' - /r')Gi?(r - r') 

r,r' 



4 — 2cOSKr — 2C0SK. 



(7) 



where Kx ~ '^{2nx + 1)/L and Hy = ky so that Gfl(r) = 
Gij(r + Ly) = —Gr{y + ix) and the charges qr obey a 
modified ncutrahty condition gr+2/ij;)mod2 = 0. 
A more convenient form of the Hamiltonian for simula- 
tion purposes is by doubling the lattice in the x direction 
to a 2L X L lattice in which the extra half is a charge 
conjugated image of the original so that 



Hr = ^VE(qr - /r)G(r - r')(gr' - /r') 



(8) 



r,r' 



where G(r) is the Green's function for & 2L x L lattice 
with periodic BC and Qr+Lx = -^r, /r+Lx = -./r @- 

To estimate the spin stiffness exponent ds, simula- 
tions were performed on a L x L lattice with eq.(^ in 
two different ways. The first is by a RT measurement 
by imposing standard periodic and antiperiodic BC cor- 
responding to /S.X = and Aj, = tt, then fitting to 



< AE^-^ >^ . This is just the procedure followed by 
all previous studies and, not surprisingly, gives essentially 
the same result = -0.76 ± 0.015 |-|J| with system 
sizes L = 4,5, 6, 7, 8, 10 and averaging over 2560 samples 
for L < 8 and 1152 for L = 10 (see Fig.(l)). This way of 
measuring a spin domain wall energy does not exploit all 
the freedom implied by eq.(||). One can find the global 
energy minimum by optimizing the BC by allowing the 
combinations {qxi — fxi) and {qyi — fyi) to vary indepen- 
dently over any integer or half integer. This corresponds 
to allowing the circulations of the phase difference and 
of Aij round the two independent loops encircling the 
torus to vary. The absolute minimum energy Eq is the 
GS energy (of a particular sample) and a spin domain 
wall is induced by — + /^^ + 1/2. The energy mini- 
mum EsD with these BC includes the energy due to the 
spin domain wall. Fitting the difference, AEf^{L) > 0, 
to eq.(|3|) yields Of^ = — 0.37± 0.015, averaging over the 
same number of samples as in the RT measurement. We 
call this a best twist (BT) measurement. This is equiva- 
lent to making a gauge transformation to all bonds in the 
direction fx = {x,y) by Aij Aij + A^/L. The energy 
E is 27r periodic in A^, E{IS.^) = E{/S.^ + 27r) and has 
a minimum at some A^ which depends on the particu- 
lar realization of disorder. The RT measurement keeps 
Jxi fixed or A^ = 0, calling the lowest energy Ep, then 
changing fxi — > fxi + 1/2 and calling the resulting lowest 
energy Eap and assuming the energy difference scales as 

. This procedure is equivalent to choosing an arbi- 
trary gauge A^{y) to compute Ep and then Eap is com- 
puted in the gauge A^ + -kS^^x/L. The original problem 
of eq.(^) is invariant under discrete gauge transforma- 
tions modulo 27r so the RT measurement is performed in 
a random gauge while the BT measurement is done in 
the gauge which minimizes the energy and depends on 
the realization of disorder. We use simulated annealing 
p2| , ^ to estimate the energy minima, which is much 
more efficient than simple quenching to T = 0. 

The chiral domain wall energy is also measured in two 
ways. Defining < AEf^ >=<]£;,„-<£;,„> | > |] 
where Em = Toim[Ep,Eap) — Er with Er the GS en- 
ergy with reflective BC gives the RT measurement for 
AE^'^ and we obtain 0^'^ = -0.37 ± 0.015. The other 
way is the BT measurement which is analogous to that 
for 9f^ when the absolute minimum energy is when the 
boundary terms in eq.(^ vanish. Since the the lowest 
energy of eq.(||) may contain a chiral but not a spin do- 
main wall, the BT condition will hold and any bound- 
ary terms must vanish. Even if, in general, there were 
boundary contributions to eq.(||), they would vanish in 
the BT condition. Thus, a BT measurement of AE^^ is 
obtained from | E^^ — E^^ \ where E^^ is the minimum 
of eq.(|) and E^'^ is the minimum of eq.(|). Fitting to 

< A£;f^(L) L^^^ yields = -0.37 + 0.010. This 
implies that Of'^ = Of^ w —0.37 to within numerical 
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accuracy, agreeing with the conjecture of Ney-Nifle and 
Hilhorst Note that the value of 6*^^^ « -0.76 does 
not satisfy the conjecture. The only difference between 
the RT and BT measurements is in Eq from eq. (o) where 
E^^ is obtained with fixed random BC and Eq hy also 
minimizing with respect to the BC. E^^ and E^^ are 
both obtained from eq.(||) and are identical because this 
is automatically a BT measurement for the special case 
of the spin glass as the boundary contributions to the en- 
ergy vanish. Note that both measurements give identical 
values for the chiral exponent 6c to within numerical un- 
certainty while the spin stiffness exponents 9f^ and 
differ by a factor of two. All 2c? results are in Fig.(l). 

Since the numerical estimates of Of^ and 9^'^ agree 
with the crucial test in 2d |^ , we can regard this as sup- 
porting our contention that we have a good definition of 
the defect energies and our numerical method is fairly 
accurate. We have done simulations on the 3dXY spin 
glass to estimate the spin stiffness exponent 9s and find 
9f'^ = +0.10 ± 0.04 with L = 2,3,4,5 (Fig.(2)). This 
is larger and more accurate than the estimate of ref. . 
The large error is due to fitting over only 3 data points. 
The negative slope of A£'f -^(i) for L = 2, 3, 4 is expected 
to become positive at larger L 1^. At present, we have 
been unable to derive the 3d analogues of eqs.(0,^, so we 
have no estimate of 9c § in 3d. 

Computations were performed at the Theoretical 
Physics Computing Facility at Brown University. JMK 
thanks A. Vallat for many discussions on XY spin glasses 
and on the importance of the CG representation when 
seeking the ground state. 
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FIG. 1. Top to bottom: L dependence of ASf^, 
ASf ^, AE^'^, A£f ^ for L = 4, 5, 6, 7, 8, 10. Solid lines 
are power law fits. 
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FIG. 2. AEf^ and ASf ^ in 3d. Solid lines are power 
law fits. 
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